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It is shown on the basis of the lowest order perturbation expansion with respect to critical 
fluctuations that the critical fluctuations give rise to an enhancement of the potential scattering 
of non- magnetic impurities. This qualitatively accounts for the enhancement of the resistivity 
due to impurities which has been observed in variety of systems near the quantum critical 
point, while the higher order processes happen to give much larger enhancement as seen from 
the Ward identity arguments. The cases with dynamical critical exponent z=2 and z=3 are 
discussed explicitly. 
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§1. Introduction 

Recently the so-called non-Fermi liquid behaviors near the quantum critical point (QCP) have 
attracted much attentions in heavy electron systems in which the tuning of the quantum parameter 
is relatively easy/"^*^^ while such phenomena have also been studied in d-band metals exhibiting 
ferromagnetic QCP.^'^"^^^ Indeed, pressures of the order of several GPa can convert the magnetic 
ground state into the paramagnetic one, or vice versa, through the magnetic QCP.^'^^^'^^'^^^ The 
anomalous behaviors of physical quantities around metallic and magnetic QCP have been success- 
fully analyzed by the so-called SCR theory (and its extentions),^'^^'^'''' except for some cases, ^'^^^ 
while the basis of SCR theory has also been established on the basis of perturbational renormaliza- 
tion group method at least up to the intermediate coupling regime. ^^'^^^ The reason why such an 
approach is successful even in the strongly correlated metals should be attributed to the validity 
of the Fermi liquid theory^'') for the description of the normal state on which the magnetic phase 
transition can be discussed. 

It has been well recognized that the impurity scattering gives rise to a drastic effect on the critical 
behaviors near the QCP.^^^ For example, the dynamical exponent z is altered considerably owing 
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to the impurity scattering: z=S in pure system for the ferromagnetic QCP is changed to z=4 
in the system with non-magnetic impurities.^^^ The effect of impurity scattering on the critical 
exponent of the temperature dependence of the resistivity has recently been discussed and attracted 
much interests.^^) However, the effect of critical fluctuations on the impurity potential has scarcely 
been discussed so far, while a general argument on the Ward identity was given in relation to the 
insulating behaviors in non-Fermi liquids subject to the impurity scattering. On the other 
hand, the enhancement of the residual resistivity po at around QCP has been reported in MnSi 
under the pressue of P 2± 15GPa^^) where the ferromagnetic state disappears, and in CeCu2Ge2 
under the pressue of P 2± 17GPa where the rapid valence change of Ce ion may occur.^'^^^ Similar 
behavior has recently been observed in UGe2 under pressures. It has also been reported that 
the resistivity of CeNi2Ge2 under P ~ 1.6GPa grows up when the temperature is decreased below 
T = 

A purpose of this paper is to discuss the effect of quantum critial fluctuations on the potential of 
non-magnetic impurity. It is shown on the basis of the perturbational treatment that the fluctua- 
tions associated with QCP can give rise to an enhancement of impurity potential which is , in the 
limit of zero-momentum transfer, in proportional to the mass enhancement factor 1/zcr, due to the 
critical fluctuations. As a result, the residual resistivity po is shown to exhibit an enhancement as 
a function of the parameter measuring quantum criticality. This consequence is quite general while 
the mass enhancement factor itself depends on the dynamical exponent, and can be compared to 
the experimental observations in a couple of systems exhibiting QCP. 

§2. Renormalization of impurity potential by critical fluctuations 
- A perturbational approach - 

In this section, we discuss the renormalization of impurity potential due to exchanging critical 
fluctuations by perturbation method. The lowest order correction of impurity potential is given 
by the Feynman diagram shown in Fig. 1. Its analytic expression for the vertex correction factor, 
corresponding to the scattering of quasiparticles from p — k/2 to p + k/2, is given as 

n' p' 

xG{p'-^,ien')G{p' + ^,ien'), (2.1) 

where A is the coupling constant between the critical fluctuation modes and the quasiparticles. The 
Green function G of quasiparticles is expressed as 

G{p,ien)= r dx^^, (2.2) 

J —oo X ICfi 

and the propagator of fluctuaion modes is expressed as 

X{q,ic0n)= r dy^^. (2.3) 
■/-oo y - IWn 
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Here, we assume that the bare quasiparticles are weh defined so that their specral weight is approx- 
imated as 

A{p,x)^z6{x-Cjj), (2.4) 

where z is the renormalization amphtude due to the effects other than the critical fluctuations; 
namely, it includes the effect of local correlation leading to heavy electron state. The spectral 
weight of fluctuation mode is given, by definition, as 

B{q, y) ^ -Imx(^, y + i5), (2.5) 
vr 

where the propagator of fluctuation modes is assumed to be parametrized as follows: 

(0) 

Xiq.u^)^ , (2.6) 

r] + A{q- Q)2 - iCgOJ 

where Xq^ is of the order of Np, the renormalized density of states (DOS) of bare quasiparticles 
at the Fermi level. The coefficient Cq in (2.6) depends on q in general and its g-dependence is 
dependent on the dynamical structure of QCP. In the case of the conventional antiferromagnetic 
(AF) fluctuations Cq is eassentially independent of q leading to the dynamical exponent z=2, while 
in the case of ferromagnetic fluctuations Cq = C/q, C being a constant, leading to z=3.^^) In the 
case of AF fluctuations where the even number of magnetic ions are contained in the unit cell and 
they are equivalent locally, ^"^^ and of uniform fluctuations which is not accompanied by conserved 
quantity^^^ such as valence fluctuations, Cq = C/max{q,£~^} in the limit T ^ 0, I being the 
mean free path of the impurity scattering. 



Fig. 1. Feynman diagram for the vertex correction of impurity potential Uq by exchanging one fluctuation mode. 



Substituting (2.2) and (2.3) into (2.1) and performing the summation with respect to n', the 
vertex correction factor for scattering potential is reduced to 



Ar,-(e + i5;^) = -E/ dyB{q,y)- 



3 



coth^ + tanhfe^^ 
2T 2T 



2T 2T 



y + i 



p+k/2-q 



iS 



(2.7) 



where the analytic continuation ie„ — ^ e + iS has been performed, and A = zA is the effective 
coupUng constant. In the hmit of T — ^ 0, using the spectral function B, (2.5), with x given by 
(2.6), the integration with respect to y in (2.7) can be easily performed obtaining 



Ar ( ^ ,^,, F(^+,^+)-F(^-,e-) 

q s-i- 

where ^± = ' ^^'^ fuctions F's are defined by 

-2E± In \E± \ + 7rsign(^±) 



(2.8) 



where 



E±ie) ^ 



Cq{e-C±) 



(2.9) 



(2.10) 



r) + Aq"^ 

Before examining the rj dependence of (2.8), let us investigate the lowest order correction for 
the selfenergy of quasiparticles due to the critical fluctuations given by (2.6). Such a selfenergy is 
given by the Feynman diagram shown in Fig. 2 which corresponds to Fig. 1 for the lowest order 
correction to the impurity potential. Analytic expression of this selfenergy is given by 



S(p,ien) = A^r^^x(^,i'^m)G(p - g,ie„ - iu>rn)- 



(2.11) 



Performing the summation with respect to m and analytic continuation ie^ — ^ e + iS, we obtain 

A 



z^ip, e + iS)= -—Y, / dyB{q, y) 

coth ^ + tanh ^ 
2T 2T 

In the limit of T — >^ 0, using (2.5) and (2.6), y-integration is easily performed obtaining 

A^ 

zS(p, e + i^) = -_ ^ x(g, 0)Fo(£:o, Cp-f), 



where the function is defined as 



Fo(^o(e),^p-g) = 



-2Eq In 1 + 7rsign(^p'_g') 
42 + 1 



where 



r] + Aq-^ 



(2.12) 



(2.13) 



(2.14) 



(2.15) 
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It is noted that the following relation holds: 



From Eqs. (2.13)-(2.15), we obtain the relation 

_ aS(p,e) _ A2 
~^ de ~ 2^ 



rj + Aq^ dEo 



(2.16) 



(2.17) 




Fig. 2. Feynman diagram for the selfenergy by exchanging one fluctuation mode. 



In the limit of forward scattering, i.e., A: — > 0, of the quasiparticles near the Fermi surface, the 
vertex correction factor (2.8) can be estimated with using (2.9) and (2.10) as follows: 

dF{E±,^±) 



\2 (J 

lim AFr Je\ ri) = — > Y(q, 0) lim 



+^Ex(9-',0)lim^-^ 



dE± 
sign(^^ 



sign(^_ 



(limfc^o^±)2 + l 



(2.18) 



The first term of (2.18) is equal to —zdTi{p,€)/d€ given by (2.17) owing to the relation (2.16). 
The second term of (2.18) vanishes, because the phase space of q satisfying the condition 
sign (^+) sign (,^_) < is restricted in a very narrow region as can be seen by geometrical con- 
sideration in the wave vector space. Indeed, its component Agy parallel to p ^ pp is restricted 
in the region |A(7||| < (k/pp)"^ because the angle between p — k/2 and p + k/2 is proportional to 
(k/pp)"^ , so that the second term of (2.18) vanishes as oc k. Therefore, we obtain 

_aS(p,e) 



l™>rr;-.(e;r/) 



(2.19) 



which implies that the renormalized impurity potential u{k) is given, in the zero momentum transfer 
limit, as 

-9S(p,e) 



u{k {)]p) 



1 



de 



u{k 0), 



(2.20) 
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where u{k) is the bare impurity potential. Namely, u is enhanced by the mass enhancement 
factor l/2;cr(p, e) = [1 — zd'E{p,e)/de)] which expresses the excess enhancement due to the critical 
fluctuations beyond the local correlations leading to the heavy electrons. This is consistent with 
the exact result obtained on the Ward identity argument by Betbeder-Matibet and Nozieres,^^^ 
who showed on the Fermi liquid formalism that the renormalized impurity potential is given as 

- 0) = ^^^Ipsf ik ^ 0), (2.21) 

where z is the renormalization amplitude including all the manybody effects and Fq the Landau 
parameter. Fermi liquid correction corresponds to that of higher order perturbation with respect to 
critical fluctuations which is beyond treatment in this paper. It is noted that renormalized impurity 
potenital depends on the momentum p of incomming quasiparticles in general, especially near the 
AF-QCP. Of the effects making z decrease, that arising from the local spin correlations should be 
cancelled by the factor (1 + Fq) in the heavy electrons as discussed by, e.g., in Ref.^°) However, 
those beyond it, such as valence or magnetic fluctuations associated with quantum criticality, can 
give rise to excess reduction of z. 

The renormalization effect of impurity potential given here is in a close relation to the conventional 
renormalization effect in the Fermi liquid theory for the conserved quantities which are expressed 
in terms of quasiparticles by the same formula as the non- interacting particles, i.e., the weight z of 
quasiparticles in the one-particle spectral weight is cancelled by the vertex correction z~^ due to 
the incoherent processes.^^'^^^ The physical reason of the enhancement of the impurity potential 
may be understood as follows: According to the expression (2.7), the vertex correction includes the 
factor B{q, y)/{y + where it is to be remembered that B{q, y) is the spectral weight of 

spin fluctuations with the wavevector q and the energy y. So, the enhancement may be related to 
that of the intermediate states of critical magnetic fluctuations associated with QCP. 

For a explicit calculation of the renormaliztion amplitude Zcr, it is convenient to rewrite (2.17) 
directly from (2.12) in the form^^) 

_5S(p,e) A2 r d?p' 



—z- 



[ ^Rex{p-p',e), (2.22) 
47r'* Jfs \Vp '\ 



de 



where Vf 's are the velocity of bare quasiparticles, and FS indicates that the integration with 
respect to p' is taken on the surface with S.f ~ £• Using the explicit form (2.6) for the propagator 
of fluctuations with Q = and the dynamical critical exponent z=3, the left hand side of (2.22) is 
calculated resulting in 



de 87r2^(uF) 



In-i^, (e = 0); 

ri 

(2.23) 

2 Aq^^+C\e\ 
13^" C\e\ ' ^^ = ')' 
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where (vp) is the averaged velocity of quasiparticles on the Fermi surface, and Qc is a cut-off 
wavenumber of the order of inverse of the lattice constant. 

Similarly, for a class of fluctuations with Q ^ and the dynamical critical exponent z=2, 
—d'S{p,e)/de given by (2.22) depends crucially on the position of p. Namely, for the momen- 
tum on the so-called "hot line" , where C^j^^q — ^Ph ~ ^ Fermi surface, essentially the 
same expression as (2.23) is obtained. Its p-dependence around the "hot line" is parameterized by 
replacing rj by rj + Aq^ in (2.23), where gm is a measure of distance from the "hot line" on the 
Fermi surface. After averaging over on the Fermi surface, one obtains 

_/_dJ:{p,e)\ 



de 



FS 



87r2A(uF) 



(e = 0); 

(2.24) 



where bi and 62 are positive constant of 0{l) depending on the details of Q and shape of the Fermi 
surface. It is remarked that the results (2.23) and (2.24) are consistent with those for the specific 
heat anomaly of SCR theroy.^^'^^^ In any case, the impurity potenial for the forward scattering is 
enhanced in proportional to l/zcr{p,e) as given by (2.20). 

§3. Effect of critical fluctuations on residual resistivity 

In order to see how this enhancement of impurity potential affects the behaviors of the resistiv- 
ity, we need to know the ^-dependence of u{k) for the scattering from p— k/2 to p -|- fc/2 near 
the Fermi surface. At first sight, the equality (2.19) holds also for general values of momentum 
transfer k because the factor [F(£^+,^+) — F(£^_, ^_)]/(^_ — ^+) in (2.8) could be approximated 
by —{dF/dE)^=^_^ for g ~ Q where x(g, 0) is diverging as I/77 near the QCP. However, this is not 
the case as shown by explicit calculation of (2.8) without using such an approximation. In this 
sense, a part of the results of Ref.^^^ should be revised as below. Indeed, the ^-dependence of u{k) 
is estimeted as follows. The expansion of [F{E+,^+) — F(£^_, ^_)]/(^_ — ^+) with respect to k in 
(2.8) is allowed so long as A; ^ \p~q\ ~ Pf, so that the relation similar to (2.18) holds. However, the 
factor limj^_^^^dF{E±,^±)/dE± in (2.18) should be replaced by dF{E±,^^_^) / dE±\^^^^^f^^_^^2^ = 
dFo{EQ, Cp-q)/^-^o|£;o=£;o(e_sjfc2), where s is a coefficient of 0(l/m) with fh being the effective mass 
of bare quasiparticles. It is because the energy of bare quasiparticle with momentum p is given by 

Ri ^py_g/2 ~ ^^"^ different from those of incoming and outgoing particles, i.e., ^^^k/2~^p-k/2^ 
where we are considering the elastic scattering. Therefore, for k ^ pp, we obtain instead of (2.19) 
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the following relation 



(3.1) 



e=e—sk^ 

Namely, the impurity potential u{k), for k <^ pp, giving the scattering on the Fermi surface (i.e., 
e=0), is renormalized as 



u{k;p) = 



1-z- 



de 



u{k;p), 



(3.2) 



e=—sk'^ 

where dT,{'p,e)/ds is given by (2.23) in the case of QCP with z=3. In the case of QCP with z=2, 
d'E{p,e)/d£ has large dependence on p as discussed above. 

If the bare impurity potential causes essentially the Born scattering, the residual resistivity po is 
enhanced by the critical fluctuations. Namely, po is given as 



po«([l + Arg.(0;??) 



x2TrNFZCiiap\u{k;p)f{l - cos9) 



(3.3) 



FS 



where Qmp is a concentration of impurity, 6 is an angle between pi k/2, and the on-shell condition 
e=^-_j_g^2=0 been used. Here it is noted that explicit dependence of renormalization amplitude 
Zcr does not appear due to cancellation between that for DOS and that for the damping rate of 
quasiparticles. We perform the calculation of (3.3) for the spherical Fermi surface with using a unit 
of wavenumber such that 2pF^/s=l, and assuming that Aq^=l. In the case of z=2, the average 
over Qia, the measure of distance from the "hot line", is approximated by J^" dqmi- • •)/ Jo" d^m- The 
results of numerical calculations of po as a function of rj are shown in Fig. 3 in an arbitrary unit 
by normalizing at 77=0. Although the logarithmic divergence in AF of (2.23) is smeared out due to 
the geometrical factor (1 — cos 9) in the case of z=3, a rather sharp cusp structure still remains in 
Po as a function of rj, like po oa 1 — Srj ioT rj 0. Here it is noted that k 2pF sm{9/2). In case of 
z=2, the result is much more smeared by the extra average process over Qj^. 



§4. Discussions 

The result obtained in §3 explains qualitatively the anomaly of po observed near the ferromagnetic 
QCP,^^'^^^ because the main anomaly arises from the factor 1/z in (2.21). In the case of AF- 
QCP, the less pronounced anomalies are expected to be observed. This is also consistent with the 
experimental fact that no strong anomaly of po is observed around AF-QCP with z=2. 

However, it needs careful consideration in the case of AF-QCP in compounds which have even 
number of magnetically equivalent ions in the primitive cell such as CeCug-ajAua;, containing mag- 
netically almost equivalent four Ce+^ ions in the primitive cell. Namely, the universality class of 
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Fig. 3. Residual resistivity po due to non-magnetic impurity as a function of inverse susceptibility rj. 



critical fluctuations of such systems belongs to that with the dynamical exponent z=3 as discussed 
previously in Ref.,^"^^ so that anomaly of po near QCP is expected to become much sharper than 
that of the conventional AF-QCP with the dynamical exponent z=2. This conclusion is consistent 
with a strong pressure dependence of po observed in CeCue,^^^ which is considered to be located 
near AF-QCP because CeCus.gAuo.i exhibits a non-Fermi liquid behavior characteristic of QCP 
although its universality class has not been identified yet. The sharp decrease of pQ as a function of 
the pressure, as shown in Fig. 2 of Ref.,^^-* does not seem to be understood as the same mechanism 
as canonical behaviors of gradual decrease of po under pressure which are observed in a series of 
heavy electron systems, such as CeInCu2,^^^ CeAls,^^) and so on. 

In the case of QCP associated with a valence transition as observed in CeCu2Ge2,^'^^ the cor- 
responding Fermi liquid effect in (2.21) may give dominant contribution and lead to much more 
pronounced enhancement of po as will be discussed elsewhere. Such a Fermi liquid correction is 
related to the higher order perturbations with respect to the critical fluctuations. Analysis of these 
higher order terms is left for future study. The Fermi liquid correction is also important for the 
enhancement of exchange potential of magnetic impruity near the ferromagnetic QCP, and leads 
to a non-trivial effect. '^^^ 

The resitivity pimp due to impurity scattering is expected to have prominent T-dependence, in 
general, arising from renormalization of impurity potential by the critical fluctuations at around 
QCP.i°) Therefore, one has to be careful when T-dependence of the observed resistivity is compared 
to existing theories.^ '''^^^ 

The enhancement of po near the ferromagneitc QCP should work to suppress the anisotropic 
supserconductivity which is expected to appear around there. One has to remember this effect 
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when dicussing the superconductivity induced by critical ferromagnetic fluctuations in real metals. 
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